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We prove that the anisotropic inflationary background of the Ackerman- Carroll- Wise model, char- 
acterized by a fixed-norm vector field, is unstable. We found the instability by explicitly solving 
the linearized equations for the most general set of perturbations around this background, and by 
noticing that the solutions diverge close to horizon crossing. This happens because one perturba- 
tion becomes a ghost at that moment. A simplified computation, with only the perturbations of 
the vector field included, shows the same instability, clarifying the origin of the problem. We then 
discuss several other models, with a particular emphasis on the case of a nonminimal coupling to the 
curvature, in which vector fields are used either to support an anisotropic expansion, or to generate 
cosmological perturbations on an isotropic background. In many cases, the mass term of the vector 
needs to have the "wrong" sign; we show that, as a consequence, the longitudinal vector mode is a 
ghost (a field with negative kinetic term, and negative energy; not simply a tachyon). We comment 
on problems that arise at the quantum level. In particular, the presence of a ghost can be a serious 
^ ■ difficulty for the UV completion that such models require in the sub-horizon regime. 

Q 

QQ ■ I. INTRODUCTION 

(— I ' Cosmological observations over the past decade have ushered in the long heralded era of precision cosmology with 
^ ^, theoretical models of the universe now being compared quantitatively to the data. In particular, Cosmic Microwave 
Background (CMB) anisotropics measurements from both orbital and sub-orbital experiments and Large Scale Struc- 
ture (LSS) observations have revealed a universe in overall agreement with the standard inflationary paradigm. This 
picture assumes an epoch preceding the standard, radiation dominated era when the scale factor increased exponen- 
^ [ tially with a nearly constant Hubble rate. This quasi deSitter expansion not only provides a solution to the classical 
'— 'i cosmological problems (homogeneity, isotropy, monopole and flatness) but also results in a nearly scale- invariant 
^ ' spectrum of super-horizon curvature perturbations which seed the structure formation process. 
, However, certain features of the full sky CMB maps recently observed by the WMAP satellite [H seem to be in 

■ conflict with the standard picture. These so called 'anomalies' include the low power in the quadrupole moment 0, 
I alignment of the lowest multipoles, also known as the 'axis-of-evil' Q, and an asymmetry in power between 

■ the northern and southern ecliptic hemispheres 6]. The statistical significance of these effects has been debated 
[ extensively in the literature. For instance, the significance of the anomalous lack of large-angle correlations, together 

psj ■ with the alignment of power has also grown in strength with the latest data 7] , with only one in 4000 realizations of 
T— I ' the concordance model in agreement with the observations. 

OO ,. Recent studies on properly masked data have shown that an anisotropic covariance matrix fits the WMAP low-i 

■ data at the 3.8a level ^ The fit was done using the full anisotropic covariance matrix and found a best-fit 
^ , anisotropy at the 10% level. The increase in significance arises from the fact that much information is contained 

; ""j ■ in the off-diagonal terms of the covariance matrix in these models llH and this is missed by simple fits to the 
' diagonal power spectrum. 

■ Such violations of statistical isotropy are considered at odds with the standard phase of early inflation. However, 
■ " " ' an albeit more plausible explanation of these anomalies arising from a systematic effect or foreground signal affecting 

the analysis is not forthcoming. This has led to a number of attempts at reconciling some of the anomalies with the 
standard inflationary picture through various modifications. 

It has been suggested in [13, [IH that the alignment of lowest multipoles could be could be due to an early anisotropic 
expansion. Indeed, an anisotropic expansion results in a nonvanishing correlation between different aem coefficients 
of the CMB multipole expansion (the formalism for dealing with the cosmological perturbations on an anisotropic 
background was also developed in [12|. fisj). To realize this with a minimal departure from the standard inflationary 
picture, the anisotropy was simply set as an initial condition at the onset of inflation. It was then shown in that 
one of the gravity waves polarization experiences a large growth during the anisotropic era (this is intimately related 
to the instability of Kasner spaces which may result in a large B signal in the CMB. Inflation however rapidly 

removes the anisotropy. The modes that leave the horizon well after the universe has isotropized were deep inside 
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^ An upper bound on the anisotropy has instead been obtained in Ref. [^. 
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the horizon while the universe was anisotropic, and one recovers a standard power spectrum at the corresponding 
scales. The signature of the earlier anisotropic stage are therefore present only on relatively large scales, which were 
comparable to the horizon when the universe was anisotropic. If the duration of inflation was too long, such scales 
get inflated well beyond the present horizon, and are irrelevant for phenomenology. Therefore, the above suggestions 
require a tuned duration of inflation. 

To avoid this tuning, one can obtain prolonged anisotropic inflationary solutions by introducing some ingredients 
that violate the premises of Wald's theorem |15| on the rapid isotropization of Bianchi universes. This has been 
realized through the addition of quadratic curvature invariants to the gravity action [16], with the use of the Kalb- 
Ramond axion [17| , or of vector fields [l^ . ^ The basic idea underlying these models is that one or more vector fields 
acquire a spatial expectation value (vev). The background solutions are homogeneous, so that the vev is only time 
dependent, and aligned along one direction; for definitcness, wc take it to coincide with the x— axis. The models differ 
from the way in which the vector field(s) acquires its vev. Most of this work we focus on one of the above models: 
namely the model by Ackerman, Carroll, and Wise (ACW) [2l[ , in which the vector acquires a fixed norm due to the 
constraint equation enforced by a lagrange multiplier. There is also a cosmological constant in the model, which is 
responsible for the accelerated expansion. The vev of the vector field controls the difference between the expansion 
rate of the x-direction and that of the two orthogonal directions. 

The original paper studied the power spectrum generated for a test field Sx in this background. It was then argued 
that this spectrum gets transferred to the spectrum of primordial density perturbations through the mechanism of 
modulated perturbations [13]. The basic idea is as follows: imagine that the effective cosmological constant in the 
model is replaced by an inflaton field, and that the decay rate of the inflaton is controlled by the vacuum expectation 
value of the test field x- Then, fluctuations Sx will result in fluctuations of the decay rate. This generates fluctuations 
of the energy density of the decay products (since the energy density of the inflaton and of the decay products scale 
differently with time). The resulting power spectrum was the starting point of the analysis of Ref. ^ cited above. The 
claim of |8|] renders the ACW model particularly interesting. However, the mechanism of modulated perturbations 
requires that the perturbations of the metric and of the energy densities which are unavoidably generated during 
inflation are subdominant; otherwise the resulting density perturbations will be a superposition of the original ones, 
generated without the test field, plus those imprinted by Sx at reheating. The original ACW work computed the 
perturbations of x on an unperturbed background. The underlying assumption is that the original perturbations in 
the metric and in the energy densities can be neglected with respect to Sx ■ In this work we show that, unfortunately, 
this assumption does not hold. On the contrary, each mode (in Fourier space) of the linearized perturbations (which 
are unavoidably generated irrespectively of the test field x) becomes nonlinear while it leaves the horizon during 
inflation. As a consequence, the background solution of [2lj is unstable, and any phenomenological study based on it 
is unreliable. 

The computation of the instability is rather complicated, due to the large number of perturbations involved. An 
important step in this computation was performed in Ref. [23|, where the linearized equations for the complete set 
of perturbations of the ACW model were presented. These equations were solved in [23] in the approximation of 
very large or very small wavelength of the perturbations. Since the Hubble expansion rates for the ACW model are 
constant, the wavelength of any mode increases exponentially with time. Therefore, at sufficiently early times during 
inflation, the wavelength is much smaller than the inverse bubble rates; this is the small wavelength regime. On the 
contrary, at sufficiently late times during inflation, the wavelength is much greater than the inverse bubble rates; this 
is the large wavelength regime. The moment at which those length become equal is known as horizon crossing We 
show here that the linearized perturbations diverge close to horizon crossing, where the equations were not solved 
in [2^ (apart from this, our computation agrees with that of [1^). This computation is presented in Section [ill It is 
just a conceptually straightforward (although technically involved) computation and explicit solution of the equations 
for the complete linearized system of perturbations for the ACW model. Each mode of the perturbations becomes 
unstable at a time close to horizon crossing (since different modes have different wavelengths, they all become 
unstable at different times). This analysis is all one needs to prove that the ACW background solution is unstable. 

To understand the reason for the instability, in Section [IIII we compute the quadratic action for the perturbations. 



^ The possibility of a late time anisotropic expansion lias instead been considered in [l^ . Although the present statistics of the observed 
Supernovae does not show any evidence for the anisotropy I20II , such studies are motivated by the large increase of data that is expected 
in the next few years. 

^ Another interesting analysis can be found in Ref. |24|| . where the late time stability of the ACW model was studied, assuming that the 

cosmological constant is replaced by a fluid. 
* The language is slightly loose, since there are two different expansion rates for the 3;— direction, and for the perpendicular y — z plane. 

However, we have in mind a small anisotropy (since a large one is forbidden by observations), so that the two rate are parametrically 

identical. 
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After the action is diagonalized, we find that the kinetic term for one of the modes goes from positive to negative 
precisely at the time t*. This has two consequences for the perturbations: 1) the equations of motion for the 
perturbations diverge at this point, as we have seen; 2) the mode is a ghost for t > t^, (we note that this second 
point alone renders the model unstable, due to the decay of the vacuum in ghost-nonghost pairs, with a rate that 
increases with the momentum of the produced particles). This last point is very interesting in the light of the findings 
of Ref. [2^. It was shown there that, in models in which the norm of the vector field is fixed as in the ACW model, it 
is possible to construct field configurations with arbitrar y la rge negative energy. Those configurations were obtained 
ignoring gravity. However, combining the result of Ref. [25| and ours, we are inclined to conclude that such models 
are generally unstable due to ghosts. 

In Sect ion HVl we present a simplified computation, with only the perturbations of the vector field included (no metric 
perturbations). This is by itself an inconsistent computation, since metric perturbations are unavoidably sourced by 
the perturbations of the vector field. However, we show it for illustrative purposes, since it is much simpler than the 
complete computations presented in the first two parts. The qualitative results of this simplified computation are in 
perfect agreement with those of the complete one (the solutions diverge close to horizon crossing, when one mode 
becomes a ghost). This can be seen with extremely simple algebra: the mode that becomes a ghost is the longitudinal 
perturbation of the vector field ("longitudinal" in the y — z plane; see Section HVl for the technical detail). Including 
the perturbations of the metric (as was done in the first two parts of this paper) simply "dresses up" the problem with 
technical complication, but does not remove the instability of the model, which is ultimately related to the vector 
field. 

The instability found for the ACW model motivated us to study also different models in which vector fields with 
broken U(l) invariance play a role during inflation. In Section |Vl we specifically focus on models in which the mass 
term is due to a nonminimal coupling of the vector to the curvature j26, 27, 28]. These models also present a ghost. 
This is not related to an anisotropic expansion (in fact, the background is isotropic for [26], in which the vector fields 
do not contribute to the expansion, and and can be made isotropic for J23] by using 3N mutually orthogonal vectors 
with equal vev), but rather to the specific sign that the mass term needs to have in these models. This is due to the 
fact that, also in this case, the ill-behaved mode is the longitudinal polarization of the vector, which is present (and, 
therefore, acquires a kinetic term) only when the mass is nonvanishing. Since this is not a trivial point, as since this 
has not yet been pointed out for the above models, we present three different proofs of this statement. In addition 
to the ghost, we have verified that also for the models [l^, the linearized equations for the perturbations diverge 
close to the horizon crossing, precisely as in the ACW model. The detailed computations will be presented elsewhere. 

In Section |Vl] we present some general discussions of these models. We do not provide quantitative computations 
here, as in the rest of the paper, but we present some comments and some worries that we believe deserve consideration. 
Particularly, we comment on some of the implications of having a ghost state in the theory. In the above models, the 
mass term provides a hard breaking of the gauge invariance. The interactions of the longitudinal vector polarization 
become strong at an energy scale parametrically set by the mass of the vector. This will occur for all the interactions 
involving the longitudinal vector, or mediated by it, renormalizing the values of the couplings between all the fields 
in the theory. Since the mass in these models is the bubble rate, or below, this may invalidate the predictions that 
rely on the initial quantization in the short wavelength regime (depending on the precise point in which the theory 
goes out of control). The most immediate solution to this problem is to provide a UV completion of the theory, in 
which the gauge symmetry is broken spontaneously from the condensation of a field that becomes dynamical above 
that mass scale. However, providing a UV completion appears to be a harder task when the longitudinal vector is a 
ghost rather than a normal field (for instance, in the higgs mechanism, the field responsible for the mass would need 
to be itself a ghost). We beheve that such issues deserve further study. 

Most of the computations presented here we summarized in our previous work [29j . 

II. INSTABILITY FROM THE LINEARIZED EQUATIONS 

In this Section we write down and solve the linearized equations for the most general set of perturbations of the 
ACW background solution. The model, and the background solution, are presented in Subsection III Al In Subsection 
III B I we introduce the perturbations; we show how to classify them in two distinct sets that are decoupled from each 
other at the linearized level, and how to construct gauge invariant variables. In Subsection III CI we write down and 
solve the linearized equations for the gauge invariant combinations. Some intermediate steps of this computation are 
put in Appendix 1X1 
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The ACW model 



The ACW model is defined by the action, 
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M, 



(1) 



where = and Vq is a constant vacuum energy (which is assumed to approximate a slow-roll inflationary 

phase). This action is a special case of the one considered in [2l[, which contains generalized kinetic terms. ^ 
The evolution of the system is governed by the field equations derived from the action ([T]): 



G lii/ 



A' 



- 2A A^ A,+g^, F' 



A (A 
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Vo 



(2) 



which are, respectively, the Einstein equations, the equation of the vector field, and the constraint forced by the 
lagrange multiplier. The constraint enforces a fixed norm of the vector. 
The background metric and vector field are 



(g^,) = diag(-l, a(^)^6(^)^M^)') 
{A^) = (0, Afpa(t) 0, 0) 



(3) 



where in the vev for the vector field {Ai) = Mp aBi we rescaled out the scale factor a and the reduced Planck mass 
Mp for algebraic convenience (since the background equations of motion are simpler when written in terms of Bi ; we 
also note that Bi is a dimensionless quantity). The last of ^ enforces a constant value for 



Mo 



(4) 



Denoting the two expansion rates Ha = a/a and i/f, — b/b, where dot is a time derivative, the nontrivial equations 
in ^ are ^ 

2Ha Hb + Hl = - Hi + 



2Hi, + 3i/2 



1 



W2 ,,2 



Hn + Hh + -ff r + Hn Hh 



1 



2 ,,2 



Vo_ 



A = H„ Hh 



I Ha 



The equations are solved by exponentially expanding scale factors a and b, with constant Hubble rates 

2Vo 1 



H^ = 



M2 6 + 7/12 + 2/14 



, Hb={l+fi^) Ha 



(5) 



(6) 



The rate Hb which characterizes the expansion of the two dimensional y ~ z plane is larger than Ha , the expansion 
rate of the anisotropic x-direction. The difference is proportional to the background expectation value of the vector 
field. The overall inflationary expansion is also supported by the vacuum energy Vb > . 



5 Specifically, the general kinetic term for the vector field considered in [lill is £ = -f^iV^ A^V - {32 (V^iA^*)^ - /Js V^A^S/^A^. Ref. 

|23| showed that the cases /3i < and /3i + & + /^s 7^ have ghosts. We have restricted our attention to the case of a standard kinetic 

term, with /3i = -fii = 1/2, ^2=0. 
® Specifically, eqs. l[5j are, respectively, the 00, 11, 22 Einstein equations, and the x— component of the vector field equation. 
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B. Classification of the perturbations and gauge invariant combinations 

We generalize to the present case the standard computation of cosmological perturbations in the case of a scalar 
field on an isotropic background (see [30| for a review) . It is convenient to exploit the symmetry of the background 
in this computation. For instance, in the isotropic case it is customary to classify the perturbations in scalar / vector 
/ tensor modes, according to how they transform under spatial rotations. The underlying reason is that modes that 
transform differently are decoupled from each other at the linearized level (namely, in the linearized equations for the 
perturbations, or, equivalently, in the quadratic expansion of the action in the perturbations). Since the modes in 
that decomposition form a complete basis for the perturbations, we could still use them here. However, this would be 
of no practical advantage, since modes belonging to different representations would now be coupled to each other. We 
can however still exploit the residual symmetry of the ACW background solution ^ in the y — z plane, and classify 
the perturbations according to how they transform with respect to rotations in this plane Modes transforming 
differently under these rotations are decoupled at the linearized level (we explicitly verify this in the equations of 
Subsection HTC)) . 

Specifically, we write the most general perturbations of the metric and of the vector field as 

/-2$ a9ix bid^B + B,) \ 

Sg^u = -2a2* abdi (d^B + B^ 

\ 62 (_2E - 2 d, dj E-d^E^- dj E,) J 

SA^ = {ao, ai, dia + ai) (7) 

where i,j = 1,2 span the isotropic coordinates. ^ The perturbations {$, x, B, ^, B, E, E, ao, ai, a} are 2d scalar 
modes: they encode 1 degree of freedom (d.o.f.) each. The perturbations {Bi, Bi, Ei, ai} are 2d vector modes. Due 
to the transversality condition {diBi = • • • = 0), they also encode 1 d.o.f. each. Notice that, contrary to the 3d 
case, there are no 2d tensor modes, since the transversality and traceless conditions eliminate all possible degrees of 
freedom. Altogether, we have 10 d.o.f. in the 2d scalar sector, and 4 d.o.f. in the 2d vector sector. These add up to 
14 d.o.f., which is the number of initial independent entries in the metric and in the vector field. 

The set of perturbations just given is redundant, since modes can be transformed into each other through infinites- 
imal coordinate transformations a;'' + f ^ . Under this change (s^l 

We decompose also the components of the infinitesimal parameter in the y — z plane in a 2d scalar plus 2d vector 
part. In this way 2d scalar (vector) modes manifestly transform into 2d scalar (vector) modes. 

To eliminate the redundancy, we can specify a gauge that completely fixes the freedom of general coordinate 
reparametrization (this was the procedure chosen in [11|, where the choice E = Y. = B = Ei^Q was made). 
Equivalently, we can construct gauge invariant combinations of the above perturbations ("gauge invariance" here 
means invariance with respect to general coordinate transformations; there is no gauge U(l) symmetry associated to 
the vector field, due to its mass term). We choose this second procedure here: each linearized equation can be written 
as "left hand side = 0"; therefore the equation must be gauge invariant. This means that the perturbations must 
enter in the linearized equations in such a way that these equations can be written in terms of the gauge invariant 
combinations only (the same is true for the quadratic action that we compute in the next Section). This provides a 
nontrivial check on our algebra. 

To find the gauge invariant combinations, we compute how each mode in ([7]) transforms under ([5]), and we then 



We chose to insert a d\ derivative in the Sqqi and 5g\i metric perturbations for algebraic convenience. For instance, we could have 
equivalently denoted 5go\ = xi- More precisely, we only study the modes of the perturbations which have a nonvanishing momentum 
component both along the x direction and in the y — z plane (in coordinate space, they have a nontrivial dependence both on x and at 
least on one between y and z). 
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arrange them into invariant combinations. We find the following set of gauge invariant modes 

$ = M„ 



* = M, 
B = 



Hb a \ a 



B- 



1 



X = 



a 



a 



1 



T. + bE 



S-a - B+-E 
a Hh \ a \ a 



1 

a 

1 

a 



a,+aMp^^-^^^^j:-bMpBidf (b + -E 



a-bMpdi {B + -E 



ao 



a 



ao-aMpBidi ( - IB + -E 



(9) 



(where the bullet denotes time derivative, and 9|n = 9| + 9|) in the 2d scalar sector, and 



B., + bE, 

a\B, + - E, 
a 

Mr, 



(10) 



in the 2d vector sector. The prefactors in front of these modes do not affect their gauge invariance, but have been 
chosen for algebraic convenience. 

Since we started from ten 2d scalars and four 2d vectors, and since there are three infinitesimal 2d scalar transfor- 
mations, and one infinitesimal 2d vector transformation, we end up with seven 2d scalar gauge invariant modes, and 
three 2d vector gauge invariant modes. One could have equivalently chosen other gauge invariant combinations of the 
initial modes. However, they can be obtained from those given here. ^ 



We perform the computations in momentum space: 



1 



(27r)^ 



-ik X 



ill) 



where 5 denotes any of the above perturbations, and where we denote with the same symbol the perturbation both in 
coordinate and in momentum space, k is the comoving momentum of the mode 5 (fc). Different modes (characterized 
each by a specific comoving momentum) are decoupled from each other at the linearized level. We denote by fc^ 
the component of the comoving momentum in the .t— direction, and by fc^ the component in the orthogonal y — z 
plane. The corresponding components of the physical momentum are p/, = ki^/a(t) and px = kx/b{t). The square 
magnitude of the comoving/physical momenta are given by k"^ — k\-\- k^ and p'^ — pj^ + Px i respectively. Finally fcri 
(pTi) denotes the component of the transverse comoving (physical) momentum in the y {i — 2) or z (i = 3) direction. 
In momentum space, the transversality conditions on the 2d vectors reads 



knVi = PTiVi = 
where Vi denotes any of the 2d vector perturbations. 



(12) 



One can verify that the present choice of variables leads to the same results as the standard computation of i30il , in the case of isotropic 
background, and of a scalar field instead of the vector. To do so, consider the gauge specified i n IllH (in practice, it amounts in identifying 
each "hatted" mode with the corresponding "non-hatted" one). It was explicitly verified in UM that, in the isotropic limit, this gauge 
leads to the same results as the longitudinal gauge. The perturbations in the longitudinal gauge can be immediately "promoted" to the 
gauge invariant modes used in [30| . This explicitly verifies that our formalism reproduces the standard one in the isotropic limit. 
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C. Solutions to the linearized equations for the perturbations 

We perturb the metric and the vector field as in eq. ([7]) . We want to solve the system of linearized equations for 
these modes. We multiply the second line of ^ by A^^ and we replace by on the right hand side. In this way 
we obtain an expression for A in terms of the vector field and the metric (this is the procedure also adopted in [3l| to 
study a similar model, and it has the advantage that we do not have to introduce explicitly the perturbation of the 
lagrange multiplier). In this way, the first and second line of ^ become 



1 „ 1 



M2 M2 



F^o. F," - A^ Vp F^P A^ A, + g^, ("J " 



9, K^^-^aA^) =0 (13) 

We then linearize these equations, and we go to momentum space as indicated in (|lip . We denote the resulting 
equations as ^ 



Eq^, : 5 [G,. ~J^T,A = , Eq^ : <5 ( — ) = (14) 



The explicit expressions are given in eqs. (|Aip in Appendix [X] The expressions pass two crucial tests: we initially 
wrote them in terms of the original perturbations ([7]) ; with some algebra, we have been able to write them solely in 
terms of the gauge invariant combinations ([9]) and pO| . As a second test, we show explicitly in Appendix [X] that the 
set of 2d scalar and 2d vector perturbations decouple in these equations. We verified that the system of 2d vectors is 
stable. For brevity, we do not report those computations here. We also show in Appendix [A] that not all the equations 
(|14p are independent (due to the perturbed Bianchi identities) . A set of independent equations is 

Eq^ : ai - ^ * = 

(15) 

Eqoo : (2 + m') ffa -I- ^ + (2 + m') (3 + 2^') Hi $ 

-2 (1 + I?) HaX-{'2 + /i^) HaB-fi Ha aa = (16) 

Eqoi: 2{l + ^l')Ha^+-^[p^J.-4fi^l + fi')Hl]x-^-2fi{l + ^i^)^ao^0 (17) 

^ 2 

Eqo- ^+(2 + /?)H„$-| + ^B-^i/„a = (18) 

Eqii : ^ + (3 + 2 ^i") hA + /^"p^ * - (2 + ^u^) i/, 4 + - (2 + fi^) (3 + 2fi^) H^] $ 

+B + (3 + 2//) HaB + + f^{5 + 2m^) Ha ao - /^Pt" = (19) 

Eqo : & + Ha& + ^ [p^ ~ 2 {l + fi^) H^] ao + fi^ ^ + fi^- ~ + A^') ^ X = 

Pt Pt Pt Pt 

(20) 

Eq- a + 3Ha& + ao+{pl-2n^Hl) a~ ^ipli' + fi^-HaB + 2Haao = Q (21) 

Pt 

where we remind that pl and px denote, respectively, the component of the physical momentum of the perturbations 
along the x direction and in the orthogonal y — z plane. We have used (fTS)) to eliminate di in favor of in all the 
other equations. Therefore, we do not need to consider this equation further in the following. 

We now solve the set of equations P^ -(PT |) . We need to rearrange them in a system that can be numerically 
integrated. There are several equivalent ways to proceed. We could either integrate all the above equations numerically, 
or we could solve some of them analytically. We choose this second option (in this way, we reduce the number of 



After we eliminated A, the linearization of the third of J2} coincides with Eqj^. 
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equations to be numerically integrated, but we generally obtain more involved expressions). We start by solving 
equations p?|) and (flS)) for x and B. First, we need to verify that this can be done. We see that this is the case, since 



a{Eqoi,EqoJ _ mMI + m')^' 



The explicit solutions are 

X 
B 



Pt 



2^2 (1 + ^2)^2 

1 p't 
2^Ji^Hl pI \l + fi- 

1 



P- 



2^PT 



Pt 



2(1+^2) 

1 

l^aPl 



Pt 



¥^0 



Pt 



ao 



Pt 



2/x (l + /i2)i?, /i 



Pt 



2nHa Pi VI + M 



9 

Pt 



(22) 



(23) 



Inserting these solutions into the remaining equations (|16p . (|20p . and (HJ), we have now a system of 4 equations 
in terms of the four unknown modes d, and do • These equations are explicitly given in appendix \X[ eqs. (|A6p 
19]) . As we show in appendix \^ by simple algebraic manipulation of these equation wc obtain the system 



( ^ ) 




(h 


do 










h 


I 4 ) 







110 

K22 K23 K24 

1 '«32 '«33 '«34 
K42 K43 K44 



(24) 



(for reasons explained in the next Section, there are no second time derivatives of do and $ in these equations). The 
entries of A4k, given in eq. (jAlOp depend on background quantities, while the four expressions at right hand side, 

given in eq. (|A11[) . are linear combinations of the unknown quantities |d, d, do, (also the coefficients of 

these linear combinations depend on background quantities). It is now straightforward to invert A^^, and to integrate 
the system numerically. 

It is easy to see that, in general, the solutions of diverge close to horizon crossing. Indeed, 



det A^. 



p^ l + /i 



Pi 



[pi - (2 + Ai') HaHb] 



(25) 



so that, when we invert this matrix, we encounter a singularity when p^ — [P' \/ Ha Hb . 

We evolve the system, starting from the adiabatic vacuum initial conditions deeply inside the horizon (eqs. (jA12|) 
and (|A16p V The two components of the physical momentum evolve as 



PL (t) 
PT (t) 



ait) 
kT 

W) 



Pto e = ptq e 



-{l+fj.'^)H^ t 



(26) 



where Plo, Pto are the values of these components when t = 0, and where we have used the background relation ^ 
between the two Hubble constants. We are free to choose the origin of the time. We set i = at the moment in which 
det = . This fixes plo = \/2 + ^2 y/Ha Hf, . The mode is then completely specified by giving the ratio pto/plo ■ 
For definiteness, we choose ^ — 0.1 (giving an O (lO~2) anisotropy), and pTo = Ha (so that the two components 
of the momentum are comparable to each other in the time range considered). We start the numerical evolution at 
t = —7/Ha, so that H/p = O (lO^'^) initially, and the modes are well inside the horizon. 

The left panel of Figure [1] shows that the modes do and $ indeed diverge when det Al^ = . The right panel gives 
the behavior close to the singularity. We give time in log units, so that t = is mapped at —log {—Ha t) = 00 . We see 
that both do and $ diverge linearly. The two modes d and 5* (not shown here) diverge logarithmically. The modes d 
and do present an identical divergence in the simplified computation of Section IIVI For that simplified analysis, we 
have also obtained the divergence analytically. It is also clear from there that the solutions diverge for generic initial 
conditions (we have also verified this in the numerical solutions of the system ([M)) ). 



9 




0.02 
0.015 
0.01 
0.005 

h 
-0.005 
-0.01 
-0.015 



Ha t . ) / agin 



8 10 

■log (- H^t) 



12 



14 



16 



FIG. 1: The right panel shows a zoom of the left panel in the region —0.6 < Hat < — 10 ^ (the time is shown in log units). 
Both modes show a 1/t divergence. 



III. GHOST FROM THE QUADRATIC EQUATION 



To gain a better understanding on the nature of the instability, we computed the quadratic action for the perturba- 
tions. We inserted the expression ^ and ([7|) in the action ([1]). We expanded at the second order in the perturbations. 
We disregarded the 2d vector modes, which are decoupled from the 2d scalars. We also eliminated the mode di 
through the constraint equation (|15p . We then Fourier transformed the perturbations as in (jlip . The final result is, 
up to boundary terms. 



<S'2dS 
'C2dS 



<fkdtah^C2dS 



2 .ap - Ha [^* * + h.c. j - (2 + Ha [<^* * + h.c.j + [B* + h.c. 



2 2 

+^i (d* * + h.c.) + ^Ha + h.c.) + (do + h.c.) - [/i^ p| + (3 + 2^^) Hi] \^\ 

-p| ($* # + h.c.) + {a* # + h.c.) - (2 + (3 + 2^l^) Hi jl-p 
+2 (1 + Ha [x* $ + h.c.) + (2 + I?) Ha (b* $ + h.c.) +^iHa (a*o^ + h.c. 



p^-2{l + ^^^) Hi 



Pl 

[pi - (3 + 2^2) ^2] |^|2 + ^ ^ + 1,.,.) 



Pi 



9 

PL 



Pl 



(27) 



The computation is conceptually straightforward, although technically involved. The resulting action passes several 
nontrivial checks. Firstly, the perturbations rearrange so that the action can be written solely in terms of gauge 
invariant combinations. Secondly, extremizing the action with respect to these modes gives precisely the linearized 
equations - (PT|) . Thirdly, the action correctly identifies the nondynamical variables of the system. The non- 
dynamical modes are the modes entering in the action (|27p without time derivatives. They do not correspond to 
propagating dynamical degrees of freedom (since the nondynamical modes can be obtained from the dynamical ones 
by their own equations of motion, that are algebraic - since they enter without time derivatives in the action). 

Whether a mode is dynamical or not depends on the initial kinetic terms. For instance, from the field strength 
Fuv = dnA^ — dyAn we immediately see that the temporal component of the vector field is nondynamical. By direct 
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inspection, we can also verify that the Ricci scalar does not lead to time derivatives for any of the (jio/x entry of the metric 
(this is most easily seen through the ADM formalism [33|). The gauge choice of [ll| {E = Ti = B = Ei = 0) preserves 
all the Sgofi modes, as well as 6Aq , namely the modes $, Xi Bi, ao- These perturbations are the nondynamical 
fields in that gauge. Therefore, the combinations x, B, Bi, cxq given in ^ and are nondynamical (since they 
are gauge invariant, and since they reduce to the nondynamical modes of the system once the gauge choice of |11| 
is made). This is explicitly verified by the quadratic action ([27)l . Wc integrate the nondynamical modes out of the 
action (j27|) . Namely, we express them as a function of the dynamical modes by using their own algebraic equations 
of motion. It is instructive to compare this with what happens in the standard case of a scalar field in an isotropic 
background. In that case, one starts from 11 perturbations (ten in the metric, and one in the scalar field). Four of 
them are removed by the gauge fixing. Four of them are nondynamical, and are integrated out (for instance, this is 
what is done in Ref. [s^] when the constraint equation (10.53) is used in going from (10.57) to (10.59)). One ends 
up with three physical propagating modes: the two polarizations of the tensor mode (gravity waves) and the scalar 
density contrast. In general it is not immediately obvious which combination of gauge invariant perturbations are 
dynamical (for instance, this is not manifest in the choice of [l^l). However, our choice of "hatted" modes makes this 
manifest. 

To be more explicit, the 2d modes entering in (P7|) can be divided in the two sets of dynamical Y = 1^", q:| and 
nondynamical N = |<i>, Xi B, doj modes. The action (P7| is formally of the type 



a., Y* Y, + (fey N* Y, + h. c.) + N* N, + {d,, Y* Y, + h. c.) + e,, Y* Y, + (/y N* Y, + h. c.) 



(28) 

where the matrices formed by the a, c, e coefficients are hermitian. The equations of motion for the nondynamical 
modes are 



SS- 



(2) 
YN 



5N* 



^ c,,N,=-h,Yj-h,Y, (29) 



These equations are precisely Eqs. p^ . (fTTI . ^5]) . and (^U)) of the linearized system that we have solved (notice 
that, indeed, they not contain any time derivative of the nondynamical modes). We solve them to express the 
nondynamical modes in terms of the dynamical ones. Inserting these expressions back into (|28p . we obtain the action 
for the dynamical modes. This action is formally of the type 



Sy^ = \ I dtd^k 



Y* K,, Y, + Y* X,, Y, + h. c. - Y* ^ Y, (30) 



where (up to boundary terms) the matrices K and OF' are hermitian, while X anti-hermitian. Extremizing this action 
we find 



5S^ 
5Y* 



= => K,, Y, + (k,, + X,, - x;) Y, + (x,, + nl) Y, = 



Yj (31) 

These two equations are precisely the linearized equations (fTO|) and (HI]), after we have inserted in them the expressions 
for the nondynamical modes. 

It is straightforwards to obtain the expressions of these matrices from the terms in the action ([27]) . However, they 
are rather involved, and not needed for the present discussion. We list instead the expression for the determinant of 
the kinetic matrix 

A ^ [2p2 + {p'+plt ~ ^Hl (1 + ^?) (2 + [2pl (1 + ^?f+pl (2 + m^) (1 + 2^^)] (32) 

One can verify that the determinant starts positive at p 3> H, and becomes zero precisely when = 
a/2 + /i^ ^JHa Hh = \/2'+^ ^/T+Ji^ Ha ■ As it is clear from the general form ([?T|) of the equations, one expect 
that the solutions diverge at this moment. This is precisely what we have found by explicitly solving these equations 
in the previous Section ! 

We also see that the determinant of the kinetic matrix becomes negative after this moment. Therefore, one of the 
physical modes of the system becomes a ghost. This by itself signals an instability of the ACW vacuum. 
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IV. SIMPLIFIED COMPUTATION 



In this Section we present a simplified stability analysis, which includes only the perturbations of the vector field. 
This analysis is not rigorous, since the metric perturbations are sourced by these perturbations already at the linearized 
level. However, as we now show, this simplified study clearly shows the origin of the instability without the involved 
algebra that is necessary to deal with the perturbations of the metric. We start from the ACW background, given in 
eqs. ([3]), (HI), and The perturbations of the vector field are as in Section HIl 

SAfj, = (ao, ai, d^a) (33) 

where we disregard the decoupled 2d vector mode a^ . 

The constraint equation enforced by the lagrange multiplier, once linearized in the perturbations, gives ai = (we 
recall, that, when the metric perturbations were included, the same constraint enforced di = /x\l/). 

The equations for the vector field (specifically, the second line of (fT3|) . in which the lagrange multiplier has been 
eliminated) then give 

- [pi - 2 (1 + ao = 

a + Haa-ao+[pl~2{l + Hi] a-Haao = (34) 

We proceed by differentiating the first equation with respect to time, and by combining it with the second equation 
so to eliminate the d term. In this way we get two first order differential equations: 

a = ^ ao 

Pt 

{l + 2^?) pi-2(l + M^) (3 + 2;.^) g,^ ^ , 
ao = p|-2(l + M^)gg i?aao-PT« (35) 

This set of equations has a singular point when p\ = 2 (l + /i^) H^ — 2 Ha Hi,. This is the divergence which 
corresponds to the one occurring at p\ = {2 + ^^) Ha H^ in the full computation of the two previous Sections. We 
set the origin of time at the singularity, so that the physical momenta are 

PL = V2(l + A*2)i?ae-^"* , pT=PToe-(i+^')^"* (36) 
Once expanded close to the singularity, the system (|35p becomes 

a w ao , ao w — ^ - Pro " (37) 

which we can combine into (we obtain the same result if we first combine the two equations (|35p into a unique 
equations for a, and we then expand that one) 

a+7+PT*a~0 (38) 

which is in turns solved be 

a « Ci Jo {pT* t) + C2 Fo {pT* t) (39) 

where Ci^2 are constants to be determined from initial conditions. While the Jq solution is constant at t = 0, the 
yj) solution has a logarithmic divergence. In principle, one may imagine arranging the initial conditions, so that the 
solution will be regular (C2 = 0) at i = 0. One would need to do so for both the real and imaginary parts of a and 
for all modes (namely, for all values of the comoving momenta) . There is however no physical reason for this tuning; 
the initial conditions are given when the mode is deeply inside the horizon, well before the equations become singular 
(there is no reason why the mode should initially "know" about the singularity that will happen close to horizon 
crossing). Moreover, as we will see, the solution is clearly divergent (C2 7^ 0) if the initial conditions are chosen in 
the adiabatic vacuum. Therefore, a has a logarithmic divergence, and ao ~ d has a linear 1/t divergence close to 
horizon crossing, which indicates that the background solution is unstable. This result is confirmed by the numerical 
solutions below. Moreover, this degree of divergence is the same as the one obtained in the complete computation of 
Subsection (|TTC|. 



12 



To find the initial conditions, and to understand the reason for the instabiUty, we compute the quadratic action for 
the perturbations. We find (up to boundary terms) 



2 



[pi+p|-2 (l + M^) H^a] l«ol 



(40) 



This actions leads to the equations (j35p . when it is extremized. The mode ao appears without any time derivatives 
in the action, therefore it can be integrated out by solving its equation of motion, which gives 



ao 



Pt 



Pi+p2,_2 (l + /i2) 
Inserting this solution back into the action we get 



a 



1 



Hi 



(41) 



(42) 



Lpi+p|,-2 (1 + Ai2) ^2 I 

In the early-time/UV limit when pj^^^ H^, this action is ghost free and a is stable. However, there is a moment later 

in the evolution, close to horizon crossing, when the longitudinal physical momentum becomes pl — (1 + /i^) Ha 
and the action vanishes. This generates the singularity in the equations (|35p . In addition, a becomes a ghost after 
the action vanishes, since the kinetic term becomes negative. 
The canonical variable of the system is 



/ pi-2(l + M^)gg 
Pi+p2,_2(l + /i2)i72 

Inserting this into the action (|42p we get (up to a boundary term) 

1 
2 



Ss 



(fkdt (|c 



(43) 



(44) 



The exact expression of lOc is rather involved, and we do not report it here. The first two terms in the early time 
expansion [H <^p) are 



_ (9 + + (3 + 2^l^yp^^ + 2 (9 + IOm^ - 4^^) pjp^ Hi f H^ 



The frequency ujc is adiabatically evolving at early times, uic <C lOc 
initial conditions for the canonical variable ar in the adiabatic vacuum 



1 



/2uv 



. dt 



. , . (45) 
O (H/p) ^ 1 . Therefore, we can set the 

(46) 



This gives the initial values of ac and dc- From (j44p and its time derivative we then obtain the initial values of a and 
d . Finally, from (|4ip we obtain the initial value of ao • Specifically, we find 



Im(ao/a). 



Pt 



1 + 



H^ 



Re (ao/a) 



M Pt Hg 
2p-^ p 



1 + 



H^ 



(47) 



The overall phase of the modes is unphysical, and we can always choose ain to be real. Since the coefRcients of (|35)) 
are real, we need to evolve this system twice, one for the real part and one for the imaginary parts of the modes. We 
did so starting from t = —7/Ha , so that H/p = O (10~^) initially. For definiteness, we take ^ = 0.1, pto = Ha (as 
in the full computation leading to Figure [T|) . We show the solutions for the real part in Figure [H The right panel is 
a close up near the singularity (which, in logarithmic units, occurs at — ln{— Hat) = +oo). The two solutions have 
been re-scaled so to make the degree of divergence manifest. The imaginary part presents an analogous divergence. 



1" The kinetic term also diverges when p = \J2 (1 + /j^) Ha (where P = yp\ +Py); this is due to the fact that oq cannot be integrated 
at this point. However, this happens after the singularity we are interested in. 
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FIG. 2: The right panel shows a zoom of the left panel in the region — 1 < Hat < —2.5 ■ 10"'^ (the time is shown in log 
units). Both modes blow up when the system (|35p . becomes singular, a oc logt and ao oc 1/t, in agreement with the analytic 
approximate solution (|39|l . Notice also that ao ^ a close to the singularity, in agreement with the analytic computation. 



V. APPEARANCE OF A GHOST IN OTHER THEORIES 



We have performed a stability analysis also of the models [l8|, |26|, |27|, |28|. The computations for these models 
are more involved than the one presented here, since they involve a larger number of perturbations. They will be 
reported elsewhere. We anticipate that these models have a ghost close to horizon crossing. In addition, for the two 
models [13, the linearized equations blow up close to horizon crossing, analogously to what happens in the ACW 
case. 

The three models [1^ [13, [11| are characterized by a nonminimal coupling of one or more vector fields to the 
curvature. The presence of the ghost is not due to an anisotropic geometry (in fact, the background evolution can be 
taken isotropic in [13], by having 3N mutually orthogonal vectors, and it is isotropic in |2Q]) but rather to the sign of 
the coupling to the curvature. The action for the vector field(s) is 



1 

12 



(48) 



where the 1/12 coefBcient allows for a slow roll of the vector in [23, [2^ (see the discussion in [23|) and the production 
of a (nearly) scale invariant spectrum of perturbations in [2^ . For simplicity, we assume that the vector field has 
no vev, and the background is pure de Sitter. Using conformal time 77, 



S= dT] (fx 



A' 



(49) 



where IvP = —2a^H^, and where the indices are risen with the inverse of the Minkowski metric 
diag(-l, +1, +1, +1). 

The fact that the longitudinal vector (|48p is a ghost in this model depends on the sign of the mass term. At first, it 
is not immediate why a mass term, and not a kinetic term, controls the nature of a mode. The main reason is that the 
longitudinal degree of freedom is present only for , since, for vanishing mass, the action is U(l) invariant, 

and the vector has only the two transverse degrees of freedom. We present three arguments to show that < leads 
to a ghost. In all the discussion below, M = |Af^| > 0. Moreover, only for this Section, = k^^ , fc^ = ki ki. 



A complete computation is presumably even more complicated for models where the kinetic term for the vector field is highly noncanon- 
ical, see for instance Ref. [33[]. 

See [331 for a similar study, based on a different model of vector fields. 

We thank A. Vainshtein for suggesting us to use the first two arguments, and for discussions. We thank A. J. ToUey for suggesting us 
to use the second argument. The third argument was presented in our previous work |29|| . 
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The most standard way to compute how many states there are in a theory, and what their nature is, is by computing 
the propagator. The quadratic lagrangian in (I49p can be cast in the form (1/2) A'' A'^ , where 



Pn 



M2 



(50) 



is the propagator. In general, the propagator needs to be diagonahzed. However, we can choose a frame in which it 
is diagonal (clearly, the number and nature of physical modes does not depend on the frame). For a positive M^, 
the pole is at = —M"^ < 0, and we can go in the rest frame, where k^^ = —k^ = (M, 0, 0, 0) . In this case, 

— {jltiu + k^j^k^/M"^^ — diag (0, —1, —1, —1), indicating that the theory has three well behaved physical particles (—1 
indicates a well behaved mode due to the signature; cf. the propagator for a scalar particle). For A'P < 0, we cannot 
go in the rest frame; however, we can choose a frame where the energy vanishes, k^ ~ k^ = (0, 0, 0, M) . In this case, 

— {'^^J.l' + kf^k^/M^) — diag(l, —1, —1, 0), indicating that one mode (the longitudinal vector) is a ghost. 

The ghost also appears in the Stuckelberg formalism. For simplicity, here M is treated has a constant (the time 
dependence of M does not modify the quadratic kinetic term, but it complicates the diagonalization) . If we redefine 



we promote the action to a gauge invariant action, with the symmetry 

1 



B„ 



(51) 



(52) 



The action ([^5]) is recovered in the unitary gauge, = 0. But one can also choose a gauge in which is transverse, 
= . In this gauge, the action of the system is 



5= / dr]d^x 



(53) 



where the field strength does not contain (p, and where the kinetic term of cj) has opposite sign to . We stress 
that the two actions ((49)) and ([53]) describe the same theory in two different gauges. We again see that the longitudinal 
component is a normal field for AP > 0, and a ghost for AP < . 

A third way to find the ghost is to decompose the vector as ~ (q;q, dia + a^) , where diai = . The action 
then separates in two decoupled pieces. 



S 



drjd'^x { - 



drjd^k 



a,- 



,'|2 



(9,aj)' 



{d.a'Y - 2d,a'd,ao - (d,af + {d.a^f + NP ao 



\OLi 



' ao + h.c.) 



(fc2 + aP^ 



ao\ 



(54) 



where prime denotes derivative wrt conformal time, and in the second line we have Fourier transformed the modes as 
in (|lip . The action splits in two decoupled parts. The first part governs the two transverse polarizations, which are 
well behaved. The second part contains only one physical mode, since ag enters without time derivatives and must 



be integrated out. The equation of motion for this field obtained from 
this back into the second part of ([54| . we obtain 



IS Oq = 



a' . Inserting 



5i 



longitudinal 



= / dr^d^k- 



fc2 M2 



/|2 



fc2 +M2 



(55) 



Again, we see that the longitudinal vector is well behaved if Al"^ > 0, and a ghost if A'P < . As a check that this 
action is correct, we can compute the equations of motion enforced by it, 

2k^aH , , 



fc2 + Af 2 



AP 



(56) 



where we have used the fact that, for Af2 = —2a?H^, one has (Af2)' — 2aHAP . This is precisely the equation of 
motion for the longitudinal vector already given in the literature, see for instance eq. (45) of (2^ (notice that, in their 
convention, R has the opposite sign wrt our convention). 
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VI. DISCUSSION 

We have shown by exphcit computations that the ACW background solution is unstable. We did so in Section |TT] 
by writing down and solving the linearized equations for the most general set of perturbations of the model. These 
equations, and their solutions, diverge at some point close to horizon crossing (this occurs for any mode in the 
Fourier decomposition; we stress that all the modes are decoupled from each other at the linearized level). The 
computations of Section |TT] prove by themselves that the background is unstable (it is possible that the divergency 
does not occur at the full nonlinear level. However, when the nonlinear interactions become important, the solution 
will be substantially different from the ACW background solution). The computations the we have presented next 
clarify the origin of the instability. Specifically, we have shown in Section IIIII that one dynamical variable of the 
system becomes a ghost at i*. The kinetic term for this mode vanishes at , which explains why the corresponding 
equation of motion, and the solutions, diverge at this point. We stress that the presence of a ghost indicates on its 
own an instability of the vacuum (see below). Our analysis is conceptually straightforward, but algebraically very 
involved. For this reason, we have also shown that our expressions pass several nontrivial checks. To gain a further 
understanding on the instability, we have also given, in Section llVi a simplified stability study, including only the 
perturbations of the vector field. The algebra involved is now considerably simpler than in the previous case. We 
verified that this system presents the same instability of the complete computation: the solutions diverge when one 
mode becomes a ghost. As we mentioned in the Introduction, the metric perturbations complicate the computation, 
but do not remove the intrinsic instability of the model, which is due to the vector field. 

In a loose sense, the ACW model is "doubly unstable" , since the linearized equations for the perturbations diverge 
- we call this "instability I" - and since it has a ghost - "instability 11" (in a generic model, either instability may 
appear without the other one). Motivated by this finding, we investigated whether other models with vector fields 
playing some role during inflation are also unstable. To our knowledge, the fist model where anisotropic inflation 
is due to a vector is the one of Ford [l^. The original proposal studied a subset of perturbations in the infinite 
wavelength limit. A general computation, performed along the line of the one presented here, shows that this model 
contains a ghost in the short wavelength regime (it has the instability II mentioned above, but not the instability I). 
The models [1^ [13, [13 are characterized by a nonminimal coupling of the vector to the curvature, with a precise 
factor. We have showed in Secton|V]why this leads to a ghost (the presence of a ghost is confirmed by the complete 
computation of the perturbations; we stress that the instability is due to the specific sign mass term, and not to the 
anisotropic expansion). The models [l^, [11] present also the instability of the type I, as defined above, although this 
is not manifest. The detailed computations will be presented elsewhere. 

An instability of the type I is a clear problem of the linearized theory. The implications of having a ghost instead 
require an understanding of the full theory. Such a study is beyond the goal of the present work. Nonetheless, we 
would like to provide some general comments. Although, contrary to what we have presented so far, these final 
considerations are not based on explicit computations, we believe that they are worth mention, and, possibly, more 
specific future investigations. 

The most obvious problem associated to a ghost is the instability of the vacuum. If a ghost is coupled to a normal 
field (and, in all the above theories, there are at least gravitational couplings), the vacuum will decay in ghost-nonghost 
excitations, with a rate which is UV divergent. This is simply due to phase space considerations. In normal situations, 
a source of any given energy can produce quanta up to some given momentum; this cuts-off the phase space available 
to the decay products. However, the quanta of a ghost field have negative energy: the higher their momentum, the 
more negative the energy is. Therefore, even a zero energy source decays into ghost-nonghost excitations conserving 
both energy and momentum. One may hope that the nonlinear interactions may somewhat arrange to cancel the 
total decay rate. However, this does not appear likely. The usual approach is to regard theories with ghosts as 
effective theories, valid only below some energy scale. Inflationary predictions heavily rely on the initial conditions; 
for instance, a slowly rolling inflaton field results in a nearly scale invariant spectrum of cosmological perturbations 
because the initial modes of the inflaton have amplitude \S(f>\ oc k~^/^ . This is due to the choice of an initial adiabatic 
vacuum in the quantized theory for the perturbations, which is performed in the deep UV regime (energy 3> H~^) [30| . 
In presence of a cut-off, the initial adiabatic vacuum cannot be imposed at arbitrarily early times, and, depending on 
the cut-off, it may not be possible to impose it at all. This casts doubts on the phenomenological predictions obtained 
for such models. 

In fact, all these theories require a cut-off which makes them invalid at high energies, irrespectively of the sign 
of the mass term. We can see this based on the behavior of massive vector fields at high energies. The models 



Several of the considerations presented here emerged from discussions with N. Kaloper, who we thank for either pointing them out, or 
for stressing their relevance. We also acknowledge very useful discussions with D.H. Lyth on these issues. 
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studied here have a gauge invariance that is broken in a hard way by the expHcit mass term for the vector. 

It is well known that, in such cases, the interactions of the longitudinal bosons violate unitarity at a scale which is 
parametrically set by M, leading to a quantum theory out of control. For the present models, M is the Hubble rate 
or below, so that the entire sub-horizon regime may be ill-defined. Although we are aware of explicit computations of 
this problem only in Minkowski spacetime, we believe that it applies also to the inflationary case, if one has unbroken 
Lorentz invariance and transitions to a locally flat frame (moreover, during inflation the mass term for the vector is 
nearly constant, and the momentum is adiabatically varying in the sub horizon regime); this issue is less transparent 
for the specific ACW model, in which the constraint enforced by A will modify the interactions of the vector; however, 
we believe that this concern is legitimate also for that model. One may then adopt the approach of ignoring the ill 
behaved longitudinal vector mode, and concentrate on the other degrees of freedom in the model. We do not believe 
that this approach is justified. At high energies, the longitudinal mode will also interact strongly with the other fields 
in the theory (this will renormalize the coupling constants). Then, depending on exactly when this happens, the 
quantum theory of the perturbations may be out of control throughout the entire short wavelength regime. If this is 
the case, all initial conditions would become unjustified, and the theory would lose its predictive power. 

We are well aware that these problems, if indeed they are present, will affect also the finding of the present work. 
However, in our opinion, the computations presented above would not lose their importance. One reason for this is 
that the instabilities that we have found have not emerged so far in the literature. Different studies have adopted 
different approximation schemes, which are very useful in the case of scalar fields, as for instance studying separately 
the short or long wavelength regimes [l^, approximating the source with fluids [35] or using the AA^ formalism 
[2^ (which allows to follow the evolution of the perturbations only in the long wavelength regime). Although such 
studies are perfectly legitimate in the regime of their validity, and can also give very useful results in the present case, 
one should be aware that these models require extra checks. We believe that our work offers some nontrivial technical 
advance in this respect. 

A more important reason is that "curing" a theory which has a hard vector mass and a ghost is more problematic 
than curing a theory with only a hard vector mass. The most immediate UV completion of a theory with a hard mass 
is through a higgs mechanism. The mass would be then due to the vev of a scalar field that becomes dynamical above 
the scale AI. In this way the theory remains under control also in the short wavelength regime, and one can apply all 
the standard computations valid for scalar fields during inflation. However, if AP needs to have the wrong sign, the 
scalar field in this UV completed theory needs to be a ghost. UV completing a theory with a condensed ghosts may 
present analogous problems to the UV completion 36] of the ghost condensation model of *37] . We hope to return to 
this issue in some future investigation. 
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APPENDIX A: EXPLICIT EXPRESSIONS FOR THE COMPUTATION OF SUBSECTION HTCl 

We write here some explicit steps in the linearized computation of the perturbations performed in Subsection III CI 



Clearly, the microscopic description of the fluids considered in |35| does not need to involve vector fields, although this has been 
considered as one of the possibilities in that work. 
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The explicit form of equations (|T4l) is 
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where the notation on the momenta is explained after eg. ljlip . and where t/j) = ^ {xi yj + Xj yi) denotes symmetriza- 
tion. Notice that it has been possible to write all the above equations in terms of gauge invariant combinations. This 
is a nontrivial check on our algebra. 

We recall that the indices i, j — 1,2 span the symmetric y — z space. Equations carrying these indices contain both 
2d scalar and 2d vector modes. However, it is easy to see that each of these equations separates in two independent 
equations, one for the 2d scalars and one for the 2d vectors: the equations in (IA1|) carrying an i index have the 
following structure 

Eqoi = kn Si + Vu , Eq^^ = kn S2 + V2i , Eq.^ = Sij S3 + kxi krj S4, + kfji V^j) , Eq^ = fc^ ^5 + (A2) 

where the expressions 5i_....5 contain the 2d scalar, and the expressions Vi^..._4 contain the 2d vector modes. Due to 
the transversality conditions of the 2d vectors. 



kn Vu = 



Vu - 



(A3) 



so that we see explicitly that the equations for the 2d scalars and the 2d vectors are indeed decoupled from each other. 
We have verified that the system of 2d vectors does not contain instabilities. For brevity, we disregard these modes 
in the following, and we concentrate only on the 2d scalars. 

We note that not all the equations in the system (IAl[) are independent. Indeed, the Bianchi identities can be 
linearized to give 



M2 \ 







(A4) 



where V^''^ denotes the covariant derivative constructed from the background metric g^i,. These can be written as 
explicit relations between the linearized equations (11411 
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We explicitly checked that the system of equations (|14p satisfies the above identities (which is a further nontrivial 
check on our algebra). Therefore, (jASp verifies that some of the linearized equations are redundant. By inserting the 
decomposition (|A2p into these identities (disregarding the 2d vector parts of the equations), we see that the i = 2,3 
identities become equal to each other. We also see that we can use the three nontrivial identities (i^ = 0, 1, 2) to 
express 52, , 53 and ^4 in terms of the other equations. Therefore, the Eqj^j and Eq^^ equations in (|Aip can be 
obtained from the other ones, and can be disregarded (obviously, one could equivalently choose to disregard some 
other equations, as long as it is possible to express them in terms of the remaining ones in (lAsp ). 

This leads us to eqs. (|15I) - (PT|) of the main text. We solve equations (fTTP and ITSl) for x and B. The other equations 
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We choose not to eliminate any more modes, but to rewrite these equations as a system of differential equations 

for XP, d, $ and do . The differential equations for do and $ are obtained by differentiating the two equations (jA6P 
and (jASp . In solving the numerical system we can replace (jA6p and (jASp with their time derivatives, provided 
that these two equations are imposed as initial conditions (see below). We are thus left with the system of eqs. 
{ (jA6P *. (|A7p . (|A8P *. (|A9P }. where (*) indicates that we take the time derivative of the equation. This system can be 
recast in the form (IMl). in terms of the coefficients 
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and 
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(All) 



As we already mentioned, we need to impose Eqs. (IA6|) and (jASp as initial condition (this allowed us to replace these 
two equations with their time derivatives). By imposing them, we obtain the initial conditions for the nondynamical 
modes in terms of the dynamical ones: 
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The initial conditions for 5*, Xf, d, ci are obtained from the early time expansion of the action (|30p . Namely, we 
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need to expand at early times the three matrices K, X, VP' entering in this action. We find 
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The next to leading order corrections to (jA13|) and (|A14[) are of O [H^/p^^ (either pL or px). Not surprisingly, the 
early time frequency is controlled by the momentum term (which is a common result in the theory of cosmological 
perturbations on a inflationary background) . What is nontrivial is the relation (|A14p between the canonical variables 
and the original modes, which, as we shall now see, is needed in order to set the initial conditions for the latter. 
Since the early time frequency varies adiabatically as 
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which are O [Ha/p) accurate. Using l|A14p we then obtain the initial conditions for the original modes as 
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The initial conditions for the nondynamical modes are then obtained by inserting these expressions into (|A12[) . These 
are the initial conditions for the numerical evolution discussed in the main text, which leads to the solutions given in 
Figure [TJ 
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